Inverse Scattering at a Fixed Energy for Long-Range 

Potentials *t 

Ricardo Weder* 

Instituto de Investigaciones en Matematicas Aplicadas y en Sistemas 
Universidad Nacional Autonoma de Mexico 
Apartado Postal 20-726, Mexico DF 01000 
weder@servidor.unam.mx 

Dimitri Yafaev 
IRMAR, Universite de Rennes I 
Campus de Beaulieu, 35042 Rennes, Cedex, France 
yafaev@univ-rennesl.fr 



Abstract 

In this paper we consider the inverse scattering problem at a fixed energy for 
the Schrodinger equation with a long-range potential in M. d , d > 3. We prove 
that the long-range part can be uniquely reconstructed from the leading forward 
singularity of the scattering amplitude at some positive energy. 



1 Introduction 



Our goal is to study the inverse scattering problem at a fixed energy for the Schrodinger 
equation with a long-range potential. We give a method for the unique reconstruction of 
the long-range part of the potential at infinity. The reconstruction of the short-range part is 
similar to the procedure suggested in ^Hj where the same problem for short-range potentials 
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was treated. To a certain extent, this paper can be considered as continuation of 13 , but 
we do not dwell here upon reconstruction of the short-range part of a long-range potential. 

For a short survey of different formulations of the inverse scattering problem see [T3] . 
Here we just mention the contributions to long-range inverse scattering. Isozaki and Kitada 
[7j consider potentials that satisfy 

\d a V{x)\ < C a (l + |a;|)- p - |Q| , p > 1/2, for all a. (1.1) 

Using stationary methods, they proved that the high-energy limit A — > oo of the scattering 
matrix S(X) determines uniquely a potential and give a method for its reconstruction. A 
similar result was obtained by Yafaev ^H] for an arbitrary p > 0. More precisely, the Fourier 
transform of the potential is reconstructed in [7] , JH] by taking the high-energy limit of the 
scattering amplitude with fixed momentum transfer. Enss and Weder |3] proved a similar 
high-energy uniqueness result and gave a reconstruction method using a time- dependent 
approach. They consider a large class of long-range potentials that are allowed to have 
singularities as well as N-body potentials. In jl] the X-ray transforms of the potentials are 
uniquely reconstructed from the action of the scattering operator on appropriate high-energy 
states. 

As mentioned above, we are interested in this paper in the inverse scattering problem at 
a fixed energy in the case of long-range potentials in M. d , d > 3. Surprisingly, this problem 
has received little attention. In |8 : it was proven that the scattering matrix at a fixed 
energy uniquely determines the asymptotics of a short-range potential on the background of 
a Coulomb potential, what actually is a different problem. We also mention paper Pj, where 
using the method of jl] the asymptotics of a long-range potential is uniquely reconstructed 
if the scattering matrix is known on some (perhaps arbitrarily small) energy interval. 

We suppose that V G C°°(M. d ), d > 3, and that, for sufficiently large \x\, 

N 

V(x)=Y,VAx) + V sr (x), (1.2) 
j'=i 

where Vj G C°°(K. d \ {0}) is a homogeneous function of order —pj, i.e., Vj(tx) = t~ Pi Vj{x) 
for all t > 0, with 1/2 < p 1 < p 2 < ■ ■ ■ < p N < 1, and V sr G C°°(R d ) is a short-range 
potential that satisfies (jl-H) for some p sr > 1. Our objective is to uniquely reconstruct 



2 



all functions Vj(x), j = 1,...,N, if only the leading forward singularity of the scattering 
amplitude is known at some fixed energy A > 0. Of course, we do not aim to reconstruct 
the whole potential, as is the case when the high-energy limit of the scattering matrix 
is known. The key issue here is that the forward singularity of the scattering amplitude 
contains all the information about the behaviour of the potential at infinity, and this allows 
us to uniquely reconstruct its asymptotic expansion at infinity. In particular, the leading 
singularity is sufficient for reconstruction of all long-range terms. Technically, we follow 
sufficiently closely our previous paper [T3], where short-range potentials that satisfy (|1-H) 
with p > 1 were considered. 

Under condition (jl.lj) the scattering matrix S(\), where A > is the energy of a quantum 
particle, is a unitary operator on L 2 (S rf_1 ). Formally, the scattering matrix can be considered 
as an integral operator, that is 



with integral kernel (the scattering amplitude) s(u,u;\). Here uo is the direction of the 
incident beam of particles and v is the direction of observation. We emphasize that our 
definition of the scattering amplitude is somewhat different from the short-range case where 
the integral kernel of S(X) and the scattering amplitude differ by the Dirac delta-function. In 
the long-range case the delta-function disappears from the integral kernel. As is well known 
PP, the scattering amplitude is C^-function away from the diagonal v — u, but its diagonal 
singularity is very wild ^3 El • 

Actually, it is more convenient (especially, in the long-range case) to consider the scat- 
tering matrix as a pseudodifferential operator. It means that 



where 11^ is the hyperplane in M. d orthogonal to u, y G Tl w is known in the physics literature 
as the impact parameter and a(y,u; A) is the right symbol of the pseudodifferential operator 
S(X). It is related to the scattering amplitude by the formula 





(1.3) 



J 



e-^a(y,uj;X)dy. 



(1.4) 
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Here and below we use freely the terminology of pseudodifferential calculus (see, e.g., [TT] . 
12 ). For example, expressions such as ()1.3|) or (|1.4|) are understood as oscillating integrals. 
Note that our definitions differ from the standard ones by the factor —k in the phase in (|1.3J) 
or (jl.4|) . The fact that the scattering matrix is well defined as a pseudodifferential operator 
was proven in [TH]. In particular, it was shown there that its symbol a(y,u; A) belongs to 
the Hormander class S° p l _ p = 5° 1 _ p (T*§ d_1 ) if p < 1 and to the to class «Si_ e e for any e > 
if p — 1. Here T*S d_1 is the cotangent bundle of the unit sphere, that is the set of points 
(y, uj) such that uj G and y G H^. The class S pl _ p fits in the standard pseudodifferential 
calculus exactly in the case p > 1/2. 

We proceed from the results of [To! ITS] where it was shown that the principal symbol 
ao(y,uj; A) of S'(A) is given by the equation 

a,{y^-\)=e- i ^' 1 ^ (1.5) 

where 

/oo 

(V(y + tu)-V(tu))dt. (1.6) 
■oo 

Roughly speaking, our approach consists of the following steps. 

1. Given S(X), we find its principal symbol ao(y, uj; A). Actually, it suffices for us to know 
the operators 

S Uo (\) = <p Uo S(\)<p Uo , (1.7) 

where uj G is an arbitrary point and ip^ is multiplication by the function (p Uo G 

C°°(S d_1 ) such that ^(w) = 1 in some, arbitrary small, neighborhood O wo of the point u . 
Then we use the fact that a (y,uj; A) coincides with the principal symbol of the pseudodif- 
ferential operator S^A) for uj G O uo . 

2. For the reconstruction of the long-range part of V, it suffices to know function (J1.6)) . 
Clearly, if V is asymptotically homogeneous of order — p, then $ is an asymptotically ho- 
mogeneous function of order — p + 1 of the variable y (except the case p — 1 when $ has a 
logarithmic behavior at infinity). Under assumption (|1.2j) the contributions to $ of different 
functions Vj can clearly be separated in (|1.6|) . Then we can directly reconstruct Vj by the 
inversion of the Radon transform (in some two-dimensional plane not passing through the 
origin). 
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Eventually, our method extends to potentials that satisfy (jl.lj) with p > 0. However, in 
the general case two new additional difficulties should be taken into account. The first is 
that the phase function in (jl.5J) is given by a more complicated formula than (jl.6|) although 
()1.6|) remains its first approximation (see ^H] > |JE] ) • The second difficulty is that for p < 1/2 
the symbol of the scattering matrix is oscillating too rapidly so that the standard pseudod- 
ifferential operator calculus cannot be applied. In this case one has to use more specific 
results for pseudodifferential operators with oscillating symbols [TEj . 

The paper is organized as follows. In Section 2, we recall different definitions of the 
wave operators in the long-range case. The scattering operator and the scattering matrix 
are also introduced there. Following ^SJ QB] > we give in Section 3 the description of leading 
forward singularity of the scattering amplitude. The classical inversion formula for the Radon 
transform is recalled in Section 4. In Section 5 we uniquely reconstruct the long-range part 
of the potential. 

2 Long-Range Scattering 

Here we recall some basic definitions of long-range scattering theory; see, e.g., ^7j, for more 
details. We consider the Schrodinger operator 

H = -A + V{x) 

with potential V(x) in the space L 2 (IR d ) where d > 2. If V is a real and bounded function, 
then the Hamiltonian H is well defined on the Sobolev class H 2 (IR d ) and is self-adjoint in 
the space L 2 (M. d ). Let us denote by Hq = —A the "free" Hamiltonian corresponding to the 
case V = 0. Under assumption the operator H has no singular continuous spectrum, 
its absolutely continuous spectrum coincides with [0, oo), and its negative spectrum consists 
of eigenvalues. 

Since the usual wave operators do not exist for p < 1, the large-time asymptotics of 
e -itH u £ Qr vec t ors ii from the absolutely continuous subspace of H is described in terms of 
the modified free evolution. There are several possibilities to construct it. For example, in 
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coordinate representation the modified free evolution Uo(t) is defined in [H] by the equation 

(U (t)u)(x) = ^ {2ity d ' 2 u{x/{2t)l (2.1) 

where 

E(x,t) = {Aty 1 \x\ 2 -t I V{sx)ds, (2.2) 

Jo 

and 

= (Fu)(0 = (27r)- rf/2 / e"^ w(x) rfx 
is the Fourier transform of u. Then the modified wave operators 

W± = s - lim e itH U (t) (2.3) 

exist and have the intertwining property HW± = W±H . Moreover, they are asymptotically 
complete, i.e., their ranges coincide with the absolutely continuous subspace of H. 

Equivalently, the modified free dynamics can be defined (see jH] and j2J) in momentum 
representation by the equation 

(FU (t)u)(O = e- l ^ 2t e-^o v ^ s)ds u(0- 

Although the operators Uo(t) and Uo(t) do not coincide, W± equals the wave operator 

W± = s- lim e itH U (t). 

Still another possibility is to define the modified free dynamics by the introduction of an 
appropriate time-independent modifier [HI EI- In this case the role of Uo(t) in the definition of 
W± is played by the operator J±e~ ltH ° where J± is a specially constructed pseudodifferential 
operator. 

Given the wave operators, the scattering operator and the scattering matrix are defined 
exactly as in the short-range case. It follows from properties of the wave operators that the 
scattering operator 

S = W* + W_ (2.4) 

commutes with H and is unitary in the space L 2 (IR d ). Let § d_1 be the unit sphere in R d , 
E+ = (0, oo) and let L 2 (R+, L 2 ^ 1 )) be the L 2 -space of functions defined on R+ with 
values in L 2 (§ d_1 ). Define the unitary operator 

T : L 2 (R d ) — > L 2 (R + ,L 2 (E> d - 1 )) 
6 



by the equation 

The spectral parameter A plays the role of the energy of a quantum particle. Then (FHqu) (A) = 
X(J r u)(X) and 

(FSu)(X) = S(X)(Fu)(X). 

The unitary operator 5(A) : L 2 (S d_1 ) — > L 2 (S rf_1 ) is known as the scattering matrix at 
energy A. 

Of course, the definition of the modified free dynamics and hence of modified wave 
operators is not unique. However, the freedom in their choice is rather limited. For example, 
in definition (j2.1j) one can add to E(x,t) an arbitrary (smooth) function which behaves as 
6±(x/(2t)) for t — > ±00. Then the wave operator ()2.3|) is replaced by W±e ld± ^ where p = 
— zV and e %6± ^ = F*e %6± ^F is multiplication by e* e± ® in the momentum representation. 
It follows that the scattering operator ()2.4j) is replaced by S = e" lS+ ' p 'Se !L ' p ' and the 
scattering matrix 5(A) is replaced by 

S{\) = e -»MVA ■) 5(A) e i M VX •). (2.5) 

In particular, if V is a sum of long-range Vi r (x) and short-range V sr (x) functions, then V 
can be replaced by Vi r (x) in (|2.2|) . In this case 

/■±oo 

^±(0 =2"' / Kr(^)^. 
>/ 

We accept below that the scattering matrix is defined in terms of wave operators (|2.3|) 
with phase function (|2.2j) . 

3 The Structure of the Scattering Matrix 

We need to know only the leading singularity of the scattering amplitude. The following 
result was essentially obtained in but it is also a consequence of more general results of 
[T5] where a complete description of all singularities was found. 
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THEOREM 3.1. Suppose that estimate (jl.lj) holds for all a. Then the scattering matrix 
S(X) is a pseudodifferential operator on the unit sphere S d_1 with the symbol 

a(y, W ; A) = e -^Hy^) (1 + %) u . A))) k = (3.1) 

where $ is function (|i.6|) and b G S^q P+1 if p < 1 and 6 G <S]~d +£ /or any e > if p = 1 . 

REMARK 3.2. Scattering matrix (j2.5|) is also a pseudodifferential operator with the sym- 
bol 

a(y, u; A) = e -*<»)- l *<*«) (1 + 6(y, «; A)), (3.2) 

where 

<%, w) = ${y, to) + 2k9+(kuj) - 2M- (kco) (3.3) 
and b belongs to the same class as the function b in Theorem 13.11 



4 Inversion of the Radon transform 

To solve our inverse scattering problem we use the two-dimensional Radon transform (see, 
e.g., [5 ). Here we briefly recall some of its properties. For v G 5~ P (1R 2 ), p > 1, the Radon 
transform, or X-ray transform, which is the same in two dimensions, is defined by the formula 

/oo 
v{ujt + y)dt 1 to G §, (uj,y) = 0. 
-oo 

It is clear that, r(y,u) = r(y, —u). The Fourier transform v of v and hence the function v 
itself can be reconstructed from its Radon transform in the following way. Let uo^ be one of 
the two unit vectors such that = 0. Hence, 



OO 



v(0 = W 1 / e-^ s r(si w € ; v)ds, £ = flfl- 1 . (4.1) 



Below we apply this method for the reconstruction of a homogeneous function V G 
C 00 (M d \ {0}) of order —p < from the integral <&(y,u) defined by formula (jl.6j) . We 
suppose that it is known for all u G S d_1 and all y G H^, y ^ 0. Actually, it suffices to know 
the function V$(w,u;) where V = V y . For an arbitrary x G M. d \ {0}, we shall find V(x). 
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Let us fix the coordinate system in such a way that the first axis is directed along x, i.e., 
x = (xi, 0, • • • , 0), and consider some two-dimensional plane A x orthogonal to x. Suppose 
that D G A x , \lu\ — 1. Differentiating (jl.6|) with respect to y±, we find that 



For y G A x such that (y, u) = 0, set 

v x (y) := d Xl V(x + y). 

Then, for all u G A x , \u\ = 1, and all y G A x , (y,uj} = 0, 

r(y,uj;v x ) = d yi $(x + y,co; V). 

Since x + y ^ 0, the function ^ G iS _p_1 (A ;c ) so that we can recover v x and, in particular, 
v x (0) = d Xl V(x) by formula (j4.1|) . Then, integrating we reconstruct 



In particular, we have proven the following proposition. 

PROPOSITION 4.1. Let V G C°°(M. d \ {0}),d > 3, be a homogeneous function of order 
-p < 0. IfV<5>{y,u) = for all u G S d_1 and ally G U w , y + 0, then V(x) = 0. 

5 Reconstruction theorem 

Reconstruction of the long-range part of the potential requires only the knowledge of the 
leading term in the asymptotics of the symbol a(y,u; A) as \y\ — > oo. The necessary result 
is formulated in Theorem 13.11 

First, we reconstruct the symbol a(y,u; A) from a family of the operators ^(A) defined 
by equation (jl.7J) . Here we use the fact that a(y,u;\) coincides with the symbol of the 
pseudodifferential operator 5 W0 (A) for u> G UQ . 

Then, by equation ()3.1|) . 




V(x) = - d s V(s,0,---,0)ds 



2ka(y,u;X) Va(y,u;\) 



zV$(y, u) + 2k(l + b(y, u; A)) _1 V6(y, u; A). 



(5.1) 
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Under condition ([1.2)1 . V6(l + b) 1 E S p where p = 2p\ > 1 if pi < 1 and p is any number 
smaller than 2 if p± = 1, whereas 

N 

W(y, cj; V) = ^ V$(y, ^) + V<%, w; Kr). (5.2) 

i=i 

Here VQ{y,uj;Vj), j = 1,---,N, are homogeneous functions of orders — pj > —1 and 
V$(-, u; V sr ) G S~ psr with — p sr < —1. Thus, given a(y,u;X), we single out in expres- 
sion ()5.1|) all homogeneous terms of orders > —1. This yields us the functions V$(y, u; Vj). 
Finally, each one of the Vj, j = 1, • ••, N, is uniquely reconstructed from V$(y, u;; Vj), as 
explained in the previous section. 

Note that the number, N, of long-range terms, as well as their order of homogeneity, 
— pj, j = 1, ■ ■ • , N, are obtained in the reconstruction process. We do not need to know them 
a priori. 

In this way we have proven our main result. 

THEOREM 5.1. Suppose that V E C°°(lR d ) ; d > 3, and that for sufficiently large \x\ 
equality JH is true. Here Vj G C°°(R d \ {0}) is a homogeneous function of order —pj with 
1/2 < pi < p 2 < ■•■ < Pn < 1j & n d V sr G C°°(M. d ) is a short-range potential that satisfies 
for some p sr > 1. Then, for an arbitrary A > 0, any family of the operators S uo (A) 
uniquely determines each Vj, j = 1, ■ ■ ■ , N. Moreover, the functions Vj can be reconstructed 
from formulae fl5.i|) . i\5. I| by the inversion of the X-ray transform. 

REMARK 5.2. Let S (A) be the scattering matrix defined by formula ()2.5|) . Considered as 
a pseudodifferential operator it has symbol ()3.2|) . According to ()3.3|) 

so that in view of ([5.2)1 the functions Vj, j = 1, . . . , N, can be reconstructed from the function 
$. Hence, we can reconstruct the long-range part of the potential from any one of the possible 
choices of scattering matrices, that correspond to different modified free dynamics. 
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